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ABSTRACT 

Thermal waves generated by either turbine rotor 
blades cutting through nonuniform combustor tempera- 
ture fields or unsteady burning could lead to thermal 
fatigue cracking In the blades. To determine the mag- 
nitude of the thermal oscillation in blades with com- 
plex shapes and material compositions, a finite 
element Galerkln formulation has been developed to 
study combustor generated thermal wave propagation In 
a model two-dimensional duct with a uniform plug flow 
profile. The reflection and transmission of the ther- 
mal waves at the entrance and exit boundaries are 
betermlned by coupling the finite element solutions at 
the entrance and exit to the eigenfunctions of an 
infinitely long adiabatic duct, Example solutions are 
presented for a variety of cases. In general, thermal 
wave propagation from an air passage into a metallic 
blade wall is small and not a problem. However, if a 
thermal barrier coating Is applied to a metallic sur- 
face under conditions of high heat transfer, a good 
impedance match is obtained and a significant portion 
of the thermal wave can pass into the blade material, 


NOMENCLATURE 

A property tensor, see Eq. < 5 > 

A t n amplitude of + going n entrance mode 

A“n amplitude of - going n entrance mode 

a n damping coefficient, Eq, (A7) 

B + n amplitude of + going exit mode 

b uniform entrance and exit duct height 

c specific heat of solid and gaseous medium at 

constant pressure 

d n propagation constant, Eq (A8) 

f frequency 


h heat transfer coefficient 

Imag Imaginary part 

i ’/T 

k thermal conductivity of the medium 

N m number of modes In expansion, Eq. (6) and (A6) 

n mode number 

n unit outward normal 

n + modified node number, Eq, (A9) 

R e real part 

t time 

U plug flow velocity 

X separation variable, Eq, (A3) 

x axial distance coordinate 

Y separation variable, Eq. (A3) 

y transverse distance coordinate 

a thermal diffusivlty, Eq. (A2) 

(3 propagation constant, Eq. (62) 

r thermal transfer coefficient, Eq. (11) 

5 element thickness at the wall 

n thermal impedance, Eq. (B9) 

e spatial (transformed) temperature 

e a entrance dimensionless temperature 

Gp exit dimensionless temperature, Eq, (7) 
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0t total dimensionless temperature. Eq, (2) 

iv separation constant, Eq, <A4) 

p density of propagation medium 

u angular velocity 

Superscript : 

~ approximate finite element prediction 

Subscripts : 

1 region 1 

2 region 2 

a entrance duct 

b exit duct 

I nodal point 

o Incident wave 

r reflected wave 

t transmitted wave 

x x direction component 

y y direction component 


INTRODUCTION 

In new aircraft combustor designs, higher outlet 
temperatures are being considered which could give 
rise to large amplitude thermal stresses In aircraft 
turbine blades and lead, perhaps, to thermal structure 
failure. Also, current earth to orbit rocket engines 
utilize hydrogen/oxygen gas generators that drive tur- 
bines used to power the oxygen and fuel pumps, These 
systems are similar to the hot section In aircraft 
turbine engines and therefore may experience similar 
dynamic thermal stresses. 

In the present paper, the sources of unsteady 
thermal waves' in the exhaust gases of a gas turbine 
engine will be discussed as well as how thermal waves 
propagate and couple to solid boundaries. Next, a 
finite element analysis Is developed which Is capable 
of evaluating thermal wave absorption In a duct with 
uniform (plug) flow and a wall composed of thermally 
absorbing materials. 


THERMAL WAVE SOURCES 

At least two mechanisms can produce thermal wave 
propagation in gas turbine passages. First, because 
of nonuniform burning In the combustor, spatial cir- 
cumferential variations In the hot gas temperature 
field exist (Norgren and Rlddlebaugh, 1985), as shown 
in Fig. 1. When the rotating turbine blades cut 
through the nonuniform circumferential temperature 
field, temporal harmonic thermal oscillations will be 
produced. For a typical turbofan engine, disturbance 
frequencies of 600 Hz or higher can be expected, 

Secondly, large temporal variations in the tem- 
perature field can exist at a fixed spatial location 


resulting from unsteady burning of the air-fuel mix- 
ture (Sofrin and Rlloff, 1976), as shown In Fig. 2, 
This temperal oscillation will effect both the turbine 
stator and rotor blades. 

The present paper Is concerned with how these 
thermal waves propagate and interact with the blade 
material . 


GEOMETRIC MODEL 

The curved turbine blades are modeled by the $ l m- 
ple rectangular geometry shown In Fig. 3(a). Next, 
the Idealized flow duct shown In Fig. 3(b) is used to 
simulate the flow In a turbine passage. The Interior 
passage of the duct Is assumed to be air moving with 
a uniform (plug) velocity profile. A metallic or 
ceramic blade represents the upper duct wall In the 
central region and the lower wall Is assumed to be a 
line of symmetry. Outside this range the entire duct 
Is assumed to be insulated. The stagnation region 
near the blade leading edge and the curvature of the 
passage have not been considered, 


METHOD OF ANALYSIS 

The thermal sources discussed In the previous 
section are assumed to generate a harmonic temperature 
field which will propagate down the duct at approxi- 
mately the flow velocity (see Appendix A, Eq. (A13>), 
As the input driving boundary condition to this prob- 
lem, a positive going thermal wave of known magnitude 
Is assumed at the entrance (x « 0) of the duct. The 
thermal wave may be plane or may have significant 
transverse temperature variations, The present paper 
will focus on the Interaction (absorption) of these 
"propagating" thermal duct waves with the wall blade 
material. High wall absorption could lead to large 
thermal stresses which could shorten the blade life 
by increasing thermal fatigue. 

In the uniform infinitely long entrance and exit 
regions with perfectly Insulated walls, the exact 
solution of the governing differential equation can 
easily be solved (see Appendix A). In contrast, exact 
analytical solutions are generally not possible In the 
central region adjacent to the turbine blade, conse- 
quently, the finite element analysis Is used to deter- 
mine the temperature fields both In the air passage 
and In the blade. 

Thermal wave reflection (usually very small) at 
the inlet and transmission (usually large) at the duct 
exit are determined by matching the finite element 
solutions In the Interior of the central region to the 
known analytical eigenfunction expansions In the uni- 
form inlet and outlet ducts. 


GOVERNING EQUATION 

The propagation of thermal disturbances will be 
governed by the convective energy equation which for 
uniform (plug) flow, neglecting viscous dissipation, 
can be written as (Bird et al,, 1960, Table 10.2-3) 



For harmonic disturbances of the form 
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For use In the finite element analysis to follow, It 
is convenient to express Eg. (3) in vector form 
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where the property tensor A is represented by 


k x 0 

0 k„ 


(4) 


(5) 


The vector product of the tensor A and a vector fol- 
low the common definition (Bird et al . , 1960, p, 730, 
Eqs. (A. 4) to ( 1 9) > , 


UNIFORM DUCT ANALYTICAL SOLUTION 

The proper termination boundary condition of the 
finite element region requires that thermal waves are 
not artificially reflected by the difference equations 
employed at the boundary nodal points (Baumelster, 
1986), The analytical solution of Eq, (3) for thermal 
wave propagation in a uniform insulated duct infinite 
in length will be employed to give the termination 
boundary condition for the finite element region. The 
analytical solution is derived in Appendix A, Eq. <A6) 
and will not be repeated here, 

Without flow, thermal waves are highly damped In 
air (see discussion In Appendix A, Eq, <A10)>. How- 
ever, when flow Is present, the thermal waves can pro- 
pagate nearly undamped in the duct like their counter- 
parts in acoustics and electromagnetics. Figure 4 
shows how only a very small Mach number is required to 
yield a unity attenuation coefficient, Thus, the pos- 
itive going wave at the entrance can be approximated 
as 

V«.,> = £ ** co S (e^ E Ua;) e , <“' l,, « (6) 

In this case, the propagation speed is equal to the 
velocity of the media, According to the usual jargon 
of acoustics and electromagnetics, the Fourier series 
terms represented by the cosine term in Eq. (6) are 
called modes. For example, for n equal to 1 , the 
transverse shape (y-dlmension) is plane and this is 
•commonly called mode 1 or simply plane wave propaga- 
tion, Mode 2 and higher are represented by the more 
complicated shapes defined by the cosine term and are 
called higher order mode propagation. In general, 
discontinuities of any kind, such as the metal blade 
in Fig, 3, will generate higher order modes because 
the uniform pattern is disturbed by the discontinuity. 


BOUNDARY CONDITIONS 

A variety of boundary conditions will be used in 
the finite element solution of Eq. (4). Each of these 
conditions will now be briefly discussed. 


Input Conditions 

The analysis assumes a given number N m of posi- 
tive going A + n waves. These waves (modes) effec- 
tively set the level of the scalar temperature field 
In the finite element region. In most of the example 
solutions to be presented, a plane wave input is 
assumed; that is, A+j is taken as unity and the rest 
of the higher modal amplitudes are assumed zero. 


Temperature Continuity 

The temperature field is continuous across an 
Interface except where heat sources are present. 
Thus, 


e a » 9 (x « 0 0 < y < b) 


(7) 


where 9 a is the modal representation of the tempera- 
ture field in the analytical inlet region given by 
Eq. (A6> of Appendix A and 9 represents the finite 
element approximation for temperature at the inter- 
face, The hat over 9 implies an approximate finite 
element numerical solution to the true solution. Sim- 
ilar equations are used at the exit. For more 
details, see Baumelster 1986, Eqs, (43) to (52). 


Heat Continuity 

Heat "flux continuity at the Interface requires 

that 


a e a 
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(inlet) 

(8) 

30b 

+k b aT 

(exit) 

(9) 


The sign change in Eqs. (8) and (9) comes directly 
from the change of the unit outward normal. 


Insulated Wal 1 

At the Insulated wall, the gradient of tempera- 
ture will be zero normal to the boundary, Thus the 
last required boundary condition is 


V§ • n « 0 


( 10 ) 


FINITE ELEMENT THEORY 


The finite element formulation of the energy 
equation was generated by using the weighted residual 
approach with the Galerkin approximation to obtain an 
Integral form of the variable property wave equation 
over the whole global domain, The continuous domain 
Is first divided into a number of discrete areas as 
shown in Fig. 3(b). This pattern as well as a symmet- 
rical diamond pattern (Baumelster 1986, Fig. 9) were 
employed. 

In the classical weighted residual manner, the 
temperature field Is curve fitted in terms of all the 
unknown modal values 9(Xf,yi>, The finite element 
aspects of converting Eq, (4) and the boundary condi- 
tions into an appropriate set of global difference 
equations can be found In text books (Sergerlfnd, 

1976) as well as Ref. 4 and for conciseness will not 
be presented herein. 
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SIMPLIFIED ANALYSIS 


Example I - Very Low Velocity Insulated Duct 


Irt Appendix B, a simple theoretical model has 
been developed to predict the thermal transfer coeffi- 
cient of a thermal disturbance In air (medium number 

1) Into a seml-inflnlte solid medium (medium number 

2) . The theory Is shown by the solid line In Fig. 

S. The simplified analysis will be useful In Inter- 
preting the more exact numerical results. The wall 
to gas temperature transfer function shown on the or- 
dinate Is a measure of the magnitude of temperature 
penetration Into the solid, The key property group 
gamma developed in Appendix B Indicates that the ther- 
mal transfer Into the blade Is proportional to 

("I 

As seen In Eq. (II), the higher the heat transfer 
coefficient, the better the energy transfer or Impe- 
dance match between the media. The transfer function 
ranges from zero (no penetration) to unity (full pene- 
tration) at the surface of the solid. The disturbance 
temperature will fall quickly in the solid material 
(see Eq. <AI0>), thus giving rise to large thermal 
stresses. 

Heat transfer coefficient on turbine blades are 
On the order of 1200 W/m 2 K (Slmoneau, 1985, Fig. 10) 
to 5700 M/m 2 K In high pressure applications (Gladden 
and Proctor, I98S, Fig. 8). With a stainless steel 
blade experiencing a driving frequency of 500 Hz, the 
Simple theory predicted a very low temperature trans- 
fer function as shown by the calculated values near 
gamma of 0.01 in Fig. 5. Consequently, thermal wave 
propagation In conventional gas turbine engines gener- 
ally have very low thermal transfer coefficients and 
does not represent a problem. 

Higher heat transfer coefficients, however, will 
lead to a better impedance match. In the SSME (Space 
Shuttle Main Engine) turbopump turbine, maximum heat 
transfer coefficients of 280 000 W/m 2 K have been pre- 
dicted (Abdul -Azt 2, 1987), Figure 5 displays trans- 
mission functions on the order of 20 percent In the 
stainless steel wall at a driving frequency of 
1000 Hz. 

Recent research (Metal Working News, 1987) at the 
NASA Lewis Research Center as well as industry has 
been concerned with ceramic coating of turbine blades 
to allow operations at higher gas temperature and 
achieve greater power plant efficiencies. In the 
schematic shown In Fig, 6, thermal barrier coatings 
(T8C) reduce the metallic blade temperature by adding 
a thin Insulator to the outer surface. The low ther- 
mal conductivity of the ceramic effectively Insulates 
the steady heat from penetrating the blade. However, 
as seen In Eq. (11), smaller values of the thermal 
conductivity k drives the value of gamma to larger 
values, A ceramic blade In the SSME environment had 
a transfer function of nearly 70 percent, as shown In 
Fig. 5. 


RESULTS OF NUMERICAL ANALYSIS 

The finite element theory will now be applied to 
some simple examples for code validation and to some 
typical hot gas turbine environments, The true value 
of the finite element theory would be In evaluating 
configurations with multiple ceramic and metallic 
layers or other geometrical complications. 


Consider the problem of an Infinitely long duct 
with Insulated walls and with a single higher order 
mode (n » 3) propagating to the right from minus 
Infinity. The amplitude of the wave is assumed unity 
at an axial position of x equals zero with a fre- 
quency of tr Hz. The ratio of Ub/a* has a value of 
50, For this case, the right going wave Is moderately 
damped as shown In Figs. 7 and 8. 

As seen In Fig. 7, the exact analytical results 
for the real and Imaginary temperatures and the finite 
element analysis are in very good agreement, The real 
and Imaginary solutions will generally display the 
spatial oscillating nature of the transformed time 
independent e(x,y) solution. Figure 8 displays simi- 
lar good agreement between the exact analytical and 
the finite element theories for the magnitude of the 
temperature oscillation. As seen in Fig. 8, the mag- 
nitude of the thermal wave decays because the thermal 
damping coefficient as predicted by the theory Is 
large Indicating that the waves energy is diffused 
Into the fluid. The magnitude of Qj ' s defined as 
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As seen In Figs. 7 and 8, the variable time does 
not appear. At a particular location, the disturbance 
temperature Sr rises and falls according to the har- 
monic forcing input conditions, defined by Eq. (2). 
Physically, the magnitude In Eq, (12> represents the 
peak value of the disturbance temperature at a given 
location at a particular time. 

In addition to these checks, the velocity, mode 
number, and fluid properties (ky> were altered to 
check the physical significants of all the terms In 
Eq. (A6) for further code validation, In all cases, 
the finite element and exact analysis were In 
agreement. 


Example 2 - Insulated Turbine Duct 

To model a typical turbojet turbine hot engine 
section, this example of the finite element technique 
considers the propagation of the a positive going 
plane thermal disturbance of 500 Hz In air moving at 
a Mach number of 0.5 In an Insulated duct without an 
absorbing duct wall. In this and all the examples to 
follow, the fluid In the main channel is assumed to 
be air at 2060 "F. Symmetrical diamond triangles were 
used to discretize the duct, The amplitude of the 
wave is assumed unity at the axial position of x 
equal to zero. The plane nature of the Impinging wave 
Is Indicated by the vertical line shown In the upper 
portion of Figs, 9 and 10, in contrast to the higher 
order mode shown in the upper portion of Fig. 7. 

Again, In this case the known analytical solution 
Eq. (A6) was In agreement with the numerical calcula- 
tions. As seen in Fig. 10, at this high Mach number 
the magnitude of the thermal wave remains at constant 
amplitude because none of Its energy Is lost to the 
Insulated wall and because the thermal damping coeffi- 
cient displayed in Fig. 4 Is near unity. 

The finite element solutions In Figs, 9 and 10 In 
the region for x/b - 0 to 1 were determined from the 
nodal values of the grid points. However, there are 
no grid points In the region for x/b greater than 1, 
as shown by the upper schematic of Fig, 9, The finite 
element values designated by the symbols for x 
greater than 1 in the lower plot of Fig, 9 or Fig, 10 
are reconstructed values from Eq, (A6) using the 
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transml tt«d mod* amplitudes predicted by the finlt* 
element program. These arc designated finite element 
solutions because the finite element analysis hasde- 
termlned the generally unknown modal reflection A“ n 
and transferred B + n coefficients. In this case, B 
replaces A In Eq. (A6) to represent the solution In 
the exit region. In the exit region, B + i equals 1 
and higher order 8 + n are nearly zero since only a 
plane wave exists. The B“ n are zero since waves 
traveling i the negative direction at the exit have 
been assumec zero. 

As seen In Fig, 9, for a frequency of 500 Hz, the 
wavelength associated with the phase nature of the 
thermal dlstjrbance is on the order of the length of 
the finite element region which was chosen to model a 
typical turbine passage. Generally 12 elements per 
oscillation would be need for accurate calculations. 
For greater accuracy more elements were employed in 
this analysis than required. 


Example 3 - 5TAINLESS STEEL DUCT HALL 

Figure 11 shows a configuration with a stainless 
steel absorber region along the upper wall, Only a 
very thin wall was required since thermal disturbances 
can only penetrate a short distance Into metallic sur- 
faces. Since there is no flow In the metallic blade, 
the wave length of the spatial oscillation will be 
proportional to a/u, as seen by the second term in 
Eq. (A10), For stainless steel, the size of the ele- 
ments In the vertical direction was required to be a 
factor of 10 smaller than In the air duct, because of 
the very small wavelength in the metallic blade, The 
axial wave propagation In the blade itself is deter- 
mined only In an approximate manner since the same 
grid spacing used In the main air duct was used in 
the metal blade, 

Figure II displays contour plots of the magnitude 
of the temperature field In the stainless steel wall, 
The maximum temperature in the blade occurs at the tip 
of the blade at x/b equals 0.1 and y/b equals 1,0. 
As shown In Fig, II, the maximum disturbance nondlmen- 
slonal wall temperature (max value of 1) Is 0,00934 
which Is very small, and quantitatively agrees with 
the analytical predictions of Fig. 5, Consequently, 
thermal stress resulting from combustor generated 
thermal waves will not be a problem In a plug flow 
configuration. 


Example 4 - Turbulent Flow Approximation 

As discussed In the Simplified Analysis section, 
the high convective heat transfer coefficient associ- 
ated with turbulent flow will produce a better energy 
transfer between the hot convective gases and the 
absorbing duct walls. An exact prediction of thermal 
disturbance propagation in turbulent flow would 
require coupling the continuity, momentum, and energy 
equations. For simplicity, the present analysis will 
now be empirically modified to roughly account for 
turbulent flow. 

Using the even grid system of Fig. 3, the trans- 
verse thermal conductivity ky In the two elements 
adjacent to the wall are now defined as 


transfer coefficients of 280 000 W/m 2 K Is used with 
an assumed driving thermal wave frequency of 1000 Hz 
produced by temporal fuel burning oscillations simi- 
lar to those shown In Fig. 2, 

In addition, to simulate the convective mixing 
of the thermal energy by the turbulent flow In the 
Central portion of the turbine duct, the ratio of con- 
ductivity in the y direction to the x direction was 
taken as 10 000. This assumption will alter the phase 
relationship In the higher order modes according to 
Eq, (A9) of Appendix A and introduces some reflections 
from the turbulent element. 

Numerical results Indicate the maximum transfer 
function of 0.2 or 20 percent of the maximum thermal 
wave temperature will enter the ceramic u:der the 
assumed conditions. This Is less than the 70 percent 
predicted by the simple theory of Appendix 6, which 
assumed a uniform bulk fluid temperature, as defined 
In Eq, (B12). However, as shown in Fig. 12, the 
finite element analysis predicts a decrease of the 
disturbance temperature field near the upper ceramic 
wall of the air duct, as shown by the disturbance tem- 
perature contour plots. The contour plot Indicates a 
bulk temperature of approximately 30 percent of the 
entrance plane wave, If the uniform mean field tem- 
perature of the simple model were reduced by 30 percent, 
the simple analytical model results and the finite 
element results would be In close agreement, Consid- 
ering the assumptions required to modify the plug flow 
theory to mode) turbulent flow, this agreement Is 
reasonable. 


CONCLUDING REMARKS 


A finite element model was developed to solve for 
the disturbance temperature field in a duct with plug 
flow. The derivations from the governing equations 
assumed that the material properties could vary with 
position resulting In a nonhomogeneous variable prop- 
erty two-dimensional thermal wave equation, This 
eliminates the necessity of finding Interface condi- 
tions between the different materials. 

Validation examples showed excellent agreement 
between finite element and exact analytical solutions. 
Numerical examples Indicated that thermal waves In 
conventional turbofan engines will not be a problem. 
However, thermal waves could possibly lead to signifi- 
cant thermal fatigue In high heat transfer regions 
such as In the SSME turbopump turbine, especially If 
thermal barrier coatings are used, 


APPENDIX A - DISTURBANCE PROPAGATION 


In the presence of uniform flow In a channel 
with insulted walls, the propagation of thermal dls 
turbances is governed by the convected energy 
Eq. (3) in the body of this report, If the fluid 
properties are assumed constant, Eq. (3) can be 
rewritten as 
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where the thermal dlffuslvlty a Is defined as 


k y - h6 (13) 

where h Is the desired heat transfer coefficient and 
the gap S represents the height of the elements. 

The heat transfer In the SSME (Space Shuttle Main 
Engine) will be used. In this case maximum heat 
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end X becomes the separation constant. For Insulated 
upper and lower walls (36/9 y - 0), the eigenfunction 
solution for Y becomes 

Y - cos rj n-1,2,3,... (A5> 

The solution for X follows the standard form 
(assume exp(mx)) to yield 


N_ 

fit 

e(x.y> » 

n* 1 


K m 




(U/2ct w - a.)x -ld_x 
le * n e 


■J7 

fit i \ 

£ A - 


(U/aa y + a n )x + ld n x 
e * " g 11 


Large Mean Ftow 

"tor large values of the mean velocity U, 

Eq. (AS) can be approximated by 

Q * A* COS 1 <w/U> X (All ) 

Most Important, the damping term has cancelled and 
thermal waves can propagate undamped In the duct like 
their analogous counter parts In acoustics and elec- 
tromagnetics. As seen In Eq. (All), the velocity of 
the thermal disturbance Is equal to the velocity of 
the medium for all modes. This contrasts with acous- 
tics, for example, for which each higher order mode 
has Its own propagation speed. 

The va'ue of the attenuation constant in Eq. <A6> 
for a plane thermal wave at a frequency of 500 Hz In 
air Is shown In Fig. 4. As seen in Fig. 4, for Mach 
numbers greater than 0.02 the attenuation factor 
approaches one (no attenuation). 


APPENDIX B - IMPEDANCE MATCHING 

The transmission of a thermal wave In a gaseous 
medium Into a solid surface will be addressed In this 
section. 
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No Flow 

Men a thermal wave Is Incident on a material 
that has properties different' from that of the medium 
in which the wave originated, part of the wave's 
energy is reflected and part' Is transmitted. This 
situation Is similar to Acoustics and Electromagnet- 
ics, Consider the case of normal Incidence of a plane 
wave In medium 1 contacting the medium 2 both of which 
are Infinite in extent at coordinate x equal to 
zero. In region I, the Incident and reflected wave 
can be expressed as 

e . e 0 e”V * e r e +l V (Bi) 

where the complex propagation constant Is equal to 
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Here, In Eq. (A6) , the A + n is the amplitude of the 
positive going waves while A“ n is the amplitude of 
negative going waves. 


No Flow 

For no convective flow (U « 0) and plane wave 
propagation to the right (n«l), Eq. (A6) reduces to 


e » a| exp 




(A10) 


13] « a| + Id] <B2) 

from Eqs . <A7> and (A8) with U set to zero. In 
region 2, only a positive going wave exists, such that 

0 « 0t e _i y <B3) 

Assuming the magnitude of the incident wave is given, 
the object of this section Is to find the magnitude of 
the reflected 0 r and transmitted 0| waves as a 
function of the thermal properties of the media. 

Continuity of the temperature at the Interface 
between the two medium requires that 

0q i* ®r “ ®t (84) 

In addition, the heat balance at the Interface 
requires that 


This is the classic solution as commonly given in text 99 

books ((Arpaci, 1986), p. 329). In contrast to elec- -k, gj 

tromagnetlc waves, thermal waves In a stationary 
medium are highly damped. For example, a thermal wave 
at a frequency of 500 Hz propagating In air at 2460 R or 

(o » 2.76 cmvsec) has a damping factor of 23.82/cm. 

Clearly, thermal waves at frequencies of 500 Hz or kj (-D]9 0 + D]9 r ) » (86) 

higher can not propagate In stationary air. 

Solving Eq. (84) and B(6> simultaneously for 0 r 

and ©t yields 
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where the thermal Impedance n Is defined as 



Thus, for Urge energy transfer the heat transfer 
coefficient should be large, Also, decrease In fre- 
quency, thermal conductivity, density or specific 
heat will improve the energy transfer. 


nj * kj|Sj 



<B9> 


The value of |J In the absence of flow is found from 
Eqs, (83) , (Al?) and <A7>, 

for no-reflectlons. nj equals n 2 , 

klP,ci « k2p 2 C2 (BIO) 

which Is independent of the driving frequency. 

for an air to metallic surface (stainless steel) 
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and in medium 2 


0 <• 0 t e e 


(613) 


where & 2 Is defined by Eq. (B2>, 

Neglecting reflection, a simple Instantaneous 
heat balance at the Interface x equal to zero 
requires that 


h,(e 0 - 0 t > - -k 2 g® 


Solving for 0* yields 

6. 


x-0 


« + k 2 |i 2 0 t 


1 


6 


o , . , h “ , . *2 


1 


1 ♦ k 


2 h, 


l + tr (1 + 1) 
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Solving for the real part of this ratio yields 


Rea 1 




^r 2 + 2 r 


+ 2 


where 


h, (Za^i/Z 

k, \ W / 
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